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A connection between the state estimation problem and the separability problem is noticed and 
exploited to find efficient numerical algorithms to solve the first one. Based on these ideas, one also 
derives a systematic method to obtain upper bounds on the maximum local fidelity when the states 
are distributed among several distant parties. 
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The estimation of specific properties of a quantum sys- 
tem subject to observation was one of the very first prob- 
lems posed in Quantum Information Theory This 
problem has no analog in the classical case, where all 
measurements commute and a sufficiently subtle observer 
(a Maxwell demon) is able to make complete tomography 
of a system without introducing any perturbations in it. 
In the quantum case, measurements are irreversible and 
invariably destroy information, and thus we have to de- 
sign and optimize our measurement strategy according 
to some figure of merit, that is, we have to decide "what 
to see" before we can extract conclusions from "what we 
see" . This figure of merit can be, for instance, the mean 
variance between our estimator and the quantity we want 
to determine 0], or the mutual information between the 
outcome of our measurement and the parameters that 
describe the possible set of states When the figure 
of merit is related to the overlap between two quantum 
states, then we are dealing with a state estimation prob- 
lem In its most usual form, this consists on estimat- 
ing the whole wave function of an unknown state when 
only a limited number of copies of this state is available. 
Initially considered as a purely theoretical problem in- 
teresting by itself, it has immediate applications in the 
calibration of quantum optical devices, where complete 
tomography is impractical, or even impossible (like when 
we deal with continuous variable systems, for example). 
Moreover, a variant of state estimation, namely, state dis- 
crimination, has recently being used in Quantum Com- 
putation to find efficient quantum algorithms to solve the 
Hidden Subgroup Problem [5|. 



be equivalent to a state estimation problem where the 
source produces orthogonal states and encodes 

them into non-orthogonal states {vf^}. A usual figure of 
merit to quantify the knowledge on provided by x is 
the so called fidelity. The idea is to prepare a state <j) x as 
a guess and then compare this state to the one that was 
originally encoded by the source just by computing the 
overlap between the two: tr^^i) (remember that we 
assume the states to be pure). If we denote by M x 

each of the Positive Operator Value Measure (POVM) el- 
ements that mathematically describe our measurement of 
the state "J/^, then the efficiency of our strategy (M x ,<f> x ) 
will then be determined by the average fidelity f: 



< / = J^PM&Mtr^i) < 1. 



(1) 



The state estimation problem consists on determining 
F, defined as the maximum fidelity among all possible 
strategies (M x ,<f> x ). 

Up to now, one has been considering global state es- 
timation. But there are many interesting situations in 
Quantum Information Theory where the encoded states 
{^ji are distributed between two or more distant par- 
ties. In this situation, the parties will have to agree on a 
strategy based on Local Operations and Classical Com- 
munication (LOCC) that allows them to extract enough 
information about the sent states so that they can make 
a good guess on the state produced by the source. Then, 
we will be interested in determining the local fidelity Fl : 



In a general state estimation scenario, a source pro- 
duces with probability pi a quantum state that is 
encoded afterwards into another quantum state 9^, to 
which we have full access. Our duty is to measure our 
given state and thus obtain a classical value x that we 
may use to make a guess on the original state vffj, which 
we will assume to be pure along this article. For abbrevi- 
ation, we will label this type of problems by (pi, ^j, \I/-). 
Games belonging to this family appear quite often in 
Quantum Information Theory. For example, we could 
think of quantum tomography as a protocol in which a 
source produces a state and encodes it into ^ = *&f N . 
Analogously, state discrimination can be easily proven to 



F L = sup£ ix p 4 tr(^(Af x ) Ai j)tr(<^* 4 ), (2) 

where the {M x )ab correspond to the POVM elements of 
a measurement accessible via LOCC. 
The main result of this paper is the following: 
Proposition 1 

For any global state estimation problem with so- 
lution F, there exists a sequence of real numbers 
pW ,F( 2 \ F( 3 \ . .., where each of the elements can be 
computed efficiently, and such that F^ > F^ > F^ > 
... > F and lim^^ F^ = F. 

This proposition can be adapted to deal with the prob- 
lem of local state estimation, and the result is 
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Proposition 2 

For any local state estimation problem with so- 
lution Fl, there exists a sequence of real numbers 
Fg , Fa , Ffp, where each of the elements can be 
computed efficiently, and such that F^ > F^ > F s 3 ^ > 
... and lim n _ 1 . 0o Fg 1 ^ — Fg > Fl, where Fg corresponds 
to the separable fidelity, i.e., the maximum attainable fi- 
delity when the POVM is separable with respect to the 
k parties. 

The structure of this paper is as follows: first, one 
will show the connection between the state estimation 
problem and the separability problem and use it to proof 
Proposition 1. An isotropic distribution of states of ar- 
bitrary dimension will be considered to illustrate the ef- 
ficiency of the resulting numerical tools. Proposition 2 
will follow naturally from Proposition 1, and the prob- 
lem of local state estimation of an arbitrary probability 
distribution of Bell states will be solved as an example. 
Finally, one will show how these ideas can be extended 
to deal with the infinite dimensional case. 

For any state estimation problem (pj, Vl/j, \E^) and any 
state estimation strategy (M x , cj> x ), the corresponding fi- 
delity / can always be expressed as 

/ = tr(pAB&AB), (3) 
with p AB = EiPi*i ® *») A AB = T, X M * ® <t>x- Note 

that, while pab is a separable quantum state fixed by 
the state estimation problem, Aab only depends on our 
measure-and-prepare strategy. It is immediate to see that 
for Aab to describe the action of a strategy over pab, 
Aab must be a separable positive semidefinite operator. 
Moreover, tr^A^s = J2 X M x = I A . Conversely, any 
separable positive semidefinite operator Aab with partial 
trace equal to the identity can be made to correspond to 
a state estimation strategy. Thus, the state estimation 
problem can be reformulated in this way: 

F = sup{tr (Aab Pab) : A A b sep, tr s A A s = I}. (4) 

To see how this approach works, consider the problem 
of estimating pure states that are distributed according 
to an isotropic probability density, i.e., we are dealing 
with (dU, U\0)(0\W, f7|0)(0|Z7+), where U is a unitary op- 
erator and dU denotes the Haar measure of the unitary 
group SU(d). Then, 

Pab — J U\0){0\tf ® U\0)(0\UUu = * + (I + V), 

(5) 

where V denotes the flip operator. Now, take Aab 
be an operator associated to a possible state estima- 
tion strategy for this problem. It is easy to see that 
/ = tr (pab Aab) = ^(pab^-ab), where A A b = J U ® 



UAabU ® WdU. Note that tr s (A AB ) = I and A A s/d G 
S, i.e., Aab corresponds to a strategy. 

The above argument shows that we can take Aab/Aa 
to be a Werner state . Werner states are a monopara- 
metrical family of states whose partial trace is propor- 
tional to the identity and whose separability regions are 
identified. Following 0], we can write Aab as Aab — 
■j^h^j [(d - t)I + (dt - 1)V], where -1 < t < 1 guaran- 
tees that Aab is positive semidefinite and t > implies 
that it is separable. A direct calculation shows that 

^{pabAab) = JTJ- ( 6 ) 

Thus, we get that F = ^-j-, recovering the results of 
Bru/3 and Macchiavello [7(. 

To solve the above problem we played with the advan- 
tage that the separability properties of Werner states are 
well known. However, most state estimation problems 
do not have the symmetries of ([5|) . How could we exploit 
the reformulation of the problem given by (j4]) in order to 
address the general case? The approach one chose makes 
use of a characterization of the set of all separable states 
due to Doherty et al. Q . This characterization is based 
on the notion of PPT symmetric extensions. 

A PPT symmetric extension A^ B to n copies of party 
A of a bipartite quantum state Aab is a positive semidef- 
inite operator A^ G B{7i t % n ®H B ) such that 

1) A^ is symmetric with respect to the interchange of 
any two copies of system A. 

2) A^g returns the state Aab when n—1 copies of system 
A are traced out. 

3) Aab remains positive under all possible partial trans- 
positions. 

We will define as the set of all bipartite operators 
Aab that admit a PPT symmetric extension to n copies 
of system A. 

If we call S to the set of all separable states, it is not 
difficult to prove that any element of S must belong to 
for any n [9(. This observation immediately leads 
us to an infinite set of necessary conditions for a state 
to be separable. Moreover, this set of conditions con- 
stitutes a hierarchy, since it follows from the definition 
of PPT symmetric extension that any state belonging to 
( 5'(™+ 1 ) na s necessarily to belong to S^ n \ What it is not 
so obvious is that this hierarchy is complete. That is, any 
quantum state belonging to for all n can be proven 
to be separable 0. 

Let dA be the dimension of 1~La- Consider now the 
following problem: 

f (n) = sup{tr(A j4B p AB ) : A AB /d A G S (n) ,tr B A AB = I}. 

(7) 

From what we have seen, it is evident that F^ > 
i?(2) > F (3) > ...> F and, moreover, lim^oo = F. 
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The advantage of this approach is that each of the bounds 
can be computed efficiently using semidefinite pro- 
gramming. 

Semidefinite programming is a branch of numerical 
analysis that is concerned with the following optimiza- 
tion problem: 

minimize c T ■ x 
subject to F{x) = F + YT l= i x i F i > °> ( 8 ) 

where {Fi}™ =0 are N x N matrices. This is known as 
the primal problem and its solution is usually denoted by 
p* . For each primal problem there is an associated dual 
problem, of the form 

d* = maximize — tr (FqZ) 
subject to Z > 0,tiFiZ = a,i = 1, ...,n, (9) 

that can also be treated a semidefinite program. It 
can be seen that d* < p* [10(. Moreover, a sufficient 
condition to assure that d* — p* is the existence of 
a strict primal feasible point, that is, a vector x such 
that F(x) > 0. Note that this condition holds in our 
case, for it is straightforward that the matrix = 
(I A /d A )® n Ols/ds > corresponds to a PPT symmet- 
ric extension to n copies of system A of the quantum 
state A AB /d A = (I A /d A ) ®l B /d B , with tv B (K AB ) = I. 

The scheme to solve dH) would then be to run algo- 
rithms that try to solve both the dual and the primal 
problem at the same time. If, after some iterations, our 
computer returns the points x. Z , we know that the so- 
lution of our problem lies somewhere in the middle, i.e., 
tr(ZFo) < p* < (F -x. Because semidefinite programs be- 
long to the P complexity class and the matrices involved 
in the calculation of each bound increase in size as 
d\d B , for fixed n, F^ can be calculated efficiently, as 
promised. We have proven Proposition 1. 

Although Proposition 1 only guarantees the asymp- 
totic convergence of the series (f ( n )) to the optimal fi- 
delity, for some specific problems we may find the value of 
F after a few iterations. Returning to the task of estimat- 
ing isotropically distributed quantum states, because the 
maximum of ^ when A AB is only demanded to be posi- 
tive semidefinite is the same as the maximum among the 
set of all separable operators, it follows that F^ = F, 

1. e., if we had performed a numerical optimization, the 
first upper bound would have coincided with the optimal 
fidelity. Also, in case the original and encoded states of 
the protocol are elements of a Hilbert space of dimension 

2, the PPT criterion is sufficient to guarantee the sep- 
arability of A AB [ll|, and thus, for all these problems, 
F = F<U. 

These ideas can be extended to obtain upper bounds on 
the solution of more complicated problems. Suppose that 
the given state is a bipartite state that is distributed 



by a source to two spatially separated parties, Alice and 
Bob. Then, Alice and Bob would have to agree on a pro- 
tocol based on local operations and classical communi- 
cation that allowed them to obtain classical information 
about the original state Vf^. 

This problem can be reformulated in a similar fashion 
as Q 

F L = sn P {ty{A ABC p A Bc) ■ ^abc LOCC}, (10) 

where p AB c = J2iPi(^i)AB ® and ^abc denotes an 
element of the set of all strategies that are accessible to 
Alice and Bob using LOCC. Unfortunately, this problem 
cannot be solved directly, for up to now no one has been 
able to characterize the set of all POVMs accessible by 
LOCC. However, we know that such a set is contained in 
the set of all separable POVMs, and thus we may conform 
with an upper bound on Fl, namely, Fs, the separable 
fidelity, which we could define as 

F s = sup{tr(A ABC p AB c) : A ABC trisep, tr c AABC = I}. 

(11) 

It is known that Fl and Fs do not coincide in general. 
For example, the separable fidelity of a uniform distribu- 
tion of the domino states [T3| would be 1, whereas Fj, can 
be proven strictly smaller than 1. However, sometimes 
F$ can be a reasonable approximation to Fl, as we will 
soon see. 

There exists a complete characterization of the set of 
all fc-separable states, also due to Doherty et al. [l3j ]. 
and also implementable using semidefinite programming 
packages. This characterization is just a generalization 
of the previous one. In this case, we would demand from 
any triseparable (normalized) state A AB q that it admits 
a PPT symmetric extension to n copies of systems A and 
B, for any n. Analogously to the global state estimation 
case, we would obtain a sequence of upper bounds Fg > 
F^ > ... > Ft,. However, this time lim, woo Fg — 
Fs > Fl- This scheme can be easily generalized to any 
number of distant parties. Proposition 2 has been proven. 

To see the usefulness of this method, consider the par- 
ticular example where Alice and Bob are distributed one 
of the four Bell states with different probabilities, and 
that each of these states encodes the very same Bell state, 
that is, if we call {"0i}i=i to the four Bell states, we would 
be dealing with the problem (pi,ipi, (ipi)AB)- Take 

V>i )2 = !(|oo)±|ii)X(oo|±(ii|), 

^3,4 — i(|01> ± |10))(<01| ± <10|). (12) 

We are going to calculate Fg 1 . The corresponding 
primal problem is 
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maximize ti(pABC^ABc) 
subject to trc(A A bc) = I, 
AabcAI a bc, a abc, a abc > 0, (13) 

with PABC = J^iPii^AB <8> (tpi)c- 

The dual of this problem (module some simplifications) 

is 

minimize tr(p) 
subject to A,B,C> 0, 
p <g>Ic — A Ta — B Tb - C Tc - pabc > 0. (14) 

From our previous discussion about semidefinite pro- 
graming, it follows that any feasible point of the dual of 
a semidefinite program provides an (in this case) upper 
bound on the solution of the primal problem. The fact 
that ij)J A = I — 2^5_i suggests that we may try to solve 
the dual problem using the following ansatz: 

A > implies that A, > 0, Vi. Analogously, p ® Ic — 
A T A — pabc > implies that 



basis. If their results are correlated, they should prepare 
the state tpi. Otherwise, they should prepare the state 
-02- It is straightforward to see that this strategy attains 
a fidelity fL=P\+P%- 

Incidentally, because the Bell states are orthogonal to 
each other, this state estimation problem is equivalent 
to the state discrimination problem (jpi, that 
is, Fl also corresponds to the maximum probability of 
determining which Bell state was sent to Alice and Bob. 

In conclusion, one has shown a systematic method to 
numerically solve the state estimation problem, and we 
have seen how this approach can be slightly modified to 
attack the local state estimation problem. Semidefinite 
programming had already been used in the resolution of 
state discrimination problems [EH], but, to the author's 
knowledge, this is the first time such a mathematical tool 
is proposed to deal with the general state estimation sce- 
nario. 

Along this paper, we have been tacitly assuming that 
we were dealing with state estimation problems in finite 
dimensions. One has ideas about how to treat the infi- 
nite dimensional case, that will be developed in future 
publications. 
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Vj ^ i. 



(16) 



What we have to minimize is the quantity / = ■ pj . 
If we fix the As, the minimum value of each pi will 
be pi = max({y : j ^ i},Pi — 4f). Therefore, / = 

J^max({-^ : j ^ i},Pi ~ -y)- Because of the symme- 
try of the problem, we can suppose that Ai > A2 are 
the two greatest As. Then, / = max(4r,pi — 4^) + 
J2j=2 max(4p)Pi — %■). This quantity can only be made 
smaller if we take A3 ; 4 = A2, and so we can assume 
/ = max(4f , Pl - 4±) + £-=2 max(^,p 4 - 4f ). 

On the other hand, for any choice of Ai, A2, it can be 
proven that / becomes strictly smaller if we take X[ — 
A' 2 = (Ai + A 2 )/2 = A. Finally, / = £•=! max(A, Pl - |). 
If we order the probabilities such that p a > pi, > p c > Pd, 
it is immediate to check that the guess A = p c provides 
the minimum value / = p a +pb- Therefore, Fl < Fjp < 
Pa +Pb- 

But such fidelity is actually attainable! All Alice 
and Bob have to do is to measure either in the com- 
putational or in the X basis each and prepare the Bell 
states with greatest probability that are compatible with 
their measurements. For example, in the case where 
Pa = Pi,Pb = P2, Alice and Bob should measure in the X 
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